Testing statistical bounds on entanglement using quantum chaos 
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Previous results indicate that while chaos can lead to substantial entropy production, thereby 
maximizing dynamical entanglement, this still falls short of maximality. Random Matrix Theory 
(RMT) modeling of composite quantum systems, investigated recently, entails an universal distribu- 
tion of the eigenvalues of the reduced density matrices. We demonstrate that these distributions are 
realized in quantized chaotic systems by using a model of two coupled and kicked tops. We derive 
an explicit statistical universal bound on entanglement, that is also valid for the case of unequal 
dimensionality of the Hilbert spaces involved, and show that this describes well the bounds observed 
using composite quantized chaotic systems such as coupled tops. 
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Recently, entanglement has been discussed extensively 
due to its crucial role in quantum computation and quan- 
tum information theory Q]. Since a quantum computer 
is a many particle system, entanglement is inevitable and 
even desirable. Entanglement is important both at the 
hardware and software levels of a quantum computer, 
as the efficiency of all proposed quantum algorithms are 
based on it, hence its characterization as a quantum re- 
source. The many particle nature of a quantum computer 
brings another phenomenon to the fore, that is chaos. 
Some studies have enquired whether chaos will help or 
hinder in the operation of a quantum computer |^ . At a 
more basic level several studies have explored the connec- 
tions between quantum entanglement and classical chaos 
1^-^ , two phenomena that are prima facie uniquely quan- 
tum and classical respectively. 

Such a connection between entanglement and chaos has 
been previously studied with the example of an TV-atom 
Jaynes-Cummings model It was found that the en- 
tanglement rate is considerably enhanced if the initial 
wave packet was placed in a chaotic region. In another 
work of similar kind, the authors have related such rates 
to classical Lyapunov exponents with the help of a cou- 
pled kicked top model . Recently, one of us studied en- 
tanglement in coupled standard maps Q and found that 
entanglement increased with coupling strength, but after 
a certain magnitude of coupling strength corresponding 
to the emergence of complete classical chaos, the entan- 
glement saturated. The saturation value depended on the 
Hilbert space dimensions and was less than its maximum 
possible value. This result implies that though there ex- 
ists a maximum kinematical limit for entanglement, dy- 
namically it is not possible to create it by using generic 
Hamiltonian evolutions on unentangled states. It should 
be emphasized that such bounds are statistical and are 
more unlikely to be violated the larger the Hilbert space 
dimension. 

Recent related work calculates the mean entangle- 



ment of pure states for the case M = N hy using a RMT 
model that allows specification of the joint probability 
distribution of the eigenvalues of the reduced density ma- 
trices (RDM) . Below we calculate the entanglement from 
a eigenvalue distribution that is valid for large M and TV. 
We show that this distribution describes well those ob- 
tained from a coupled kicked top model. There is also 
some early work that calculates the subsystem entropy 
for random pure states 0. Apart from RMT simulations, 
we deal with an actual quantum mechanical system and 
relate these results to the presence of classical chaos, thus 
our results show in what context results such as in 
can be expected to be universal. 

The previous studies on entanglement, in the context 
of chaos, were based on pure states of bipartite systems, 
where the von Neumann entropy of the RDM is a nat- 
ural measure of quantum entanglement. We will also 
initially consider pure states and point out in the end 
that as a simple corollary we can estimate the entangle- 
ment of formation of any density matrix as well. Sup- 
pose that the state space of a bipartite quantum system 
isTi. — Hi® H2, where dimHi = N < dimH2 = M, and 
dimTi — d = NM. If p = J2iP'i\'f'i){4'i\ is an ensemble 
representation of an arbitrary state in 7i, the entangle- 
ment of formation is found by minimizing X^i Pi-^(l'^j)) 
over all possible ensemble realizations. Here E is the von 
Neumann entropy of the RDM of the state \(j>i) belonging 
to the ensemble, i.e., its entanglement. For pure states 
\^) there is only one unique term in the ensemble rep- 
resentation and the entanglement of formation is simply 
the von Neumann entropy of the RDM. 

The two RDMs of the bipartite state \tp) are pi — 
Tr2(|i/')(V'I) and = TTi{\ip) {ip\) . The Schmidt decom- 
position of IV') is the optimal representation in terms of 
a product basis and is given by 
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where < Ai < 1 are the (nonzero) eigenvalues of either 
RDMs and the vectors are the corresponding eigenvec- 
tors. The von Neumann entropy Sy is the entanglement 
E{\^)) given by 



U 



exp 



N 



TrKpilnpO = -2^A,ln(A,),; 1 = 1,2. (2) 



Under an arbitrary unitary evolution p or \^) evolves 
into states with changed entanglement. Quantum chaotic 
evolutions eventually create large entanglement which 
fluctuates around the value ln(7A^). The factor 7 de- 
pends only on the ratio M/N and tends to unity (maxi- 
mal entanglement) as M — > cx). Such evolutions lead to 
universal properties of the RDMs, which is also shared 
by the RDMs of stationary states. This follows from the 
universal near normal distributions of the complete pure 
state components in any generic bases. The distribu- 
tion of the eigenvalues of RDMs {A^} also follows from 
RMT results for correlation matrices recently used in the 
analysis of data from financial time series . Many im- 
portant universal spectral fluctuation properties of quan- 
tum chaotic systems have been modeled and explained by 
RMT. We extend this success to the RDMs of compos- 
ite systems and consequently a description of quantum 
entanglement in strongly interacting systems. 

As the Hilbert space dimension and chaos have roles 
in this bound for entanglement, coupled large spins are 
attractive models. A coupled kicked tops model has al- 
ready been used in this context [Q, we generalize it here 
to include the case of unequal spins and symmetry break- 
ing terms. The Hamiltonian of the coupled top system 
used is : 

H{t) = \jy, + ^^{J^,+a,f '^(^-") 
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The Jy^ terms describe free precession of each top and the 
remaining terms are due to periodic (5-function kicks. The 
first two such terms are torsion about z— axis and the fi- 
nal term describes the spin-spin coupling. When either of 
the constants, ai or a2, is not zero the parity symmetry 
RH{t)R~^ = H{t), where R = exp(i7rJj^J (g) exp{iTrJy^), 
is broken. The dimensionality of the Hilbert spaces are 
N = 2ji-|-l and M = 2^2 + 1. The unitary time evolution 
operator corresponding to this Hamiltonian is given by : 

Ut = (c/i ® U2) = [ {uiuD ® (uiu^) ] UI2, (4) 

where the different terms are given by 
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and i = 1,2. There exists an antiunitary general- 
ized time reversal symmetry, [exp(i7r J^j ) exp^in Jy-^ /2)](E) 
[cxp{iTTJx2 ) exp^inJy^ where K is complex conjuga- 
tion operator, from which we can expect the applicability 
of results concerning the Gaussian orthogonal ensemble 
(GOE). We note that for the parameter values consid- 
ered in this Letter, the nearest neighbor spacing distri- 
bution (NNSD) of the eigenangles of Ut is Wigner dis- 
tributed, which is typical of quantized chaotic systems 
with time reversal symmetry. Entanglement production 
of time evolving states under Ut have been studied for 
two different initial states. (1) The initial state is a prod- 
uct of directed angular momentum states as given in Rcf. 
[Q, placed in the chaotic sea of phase space. This is a 
completely unentangled state. (2) The initial state is 
maximally entangled and is given by: 



(mi,TO2|V'(0)) = ^S,nim2- 



(6) 



These initial states are evolved under Ut, and the results 
are displayed in Fig. |l|. Here the coupling strength is very 
large compared to the value taken in 0| , as our goal is to 
study entanglement saturation and strong coupling will 
help us achieve entanglement saturation within a short 
time. In the first case, initially both the von Neumann 
entropy and the linearized entropy (Sr = 1 — Tri(p^)), 
are zero, but with time evolution both entropies start 
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FIG. 1. Entanglement saturation of a completely unentan- 
gled initial state (solid line) and a maximally entangled initial 
state (dotted line) under time evolution operator Ut- Here 
fc = 3, e = 0.1 and the phases 0:1=02 = 0.47. Inset shows 
similar behaviour of linear entropy. 
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increasing and get saturated, apart from small fluctua- 
tions, at values less than their maximum possible values. 

For the von Neumann entropy the saturation value is 
~ ln(0.67V) and for the linear entropy it is approximately 
1 — 2/N , where N is the dimension of each subsystems. 
This is the dynamical bound for entanglement of a sys- 
tem consists of two equal dimensional subsystems, while 
the maximum kinematical limits are In TV and 1 — 
respectively. The saturation value of von Neumann en- 
tropy of this time evolved state is same as that obtained 
for stationary states of completely chaotic coupled stan- 
dard maps [^. In the second case, the initial state is 
maximally entangled and time evolution forces this state 
to partially disentangle till the entropy reaches the above 
mentioned values. 

This study shows that the saturation of entanglement 
is a universal phenomenon, it depends only on the Hilbert 
space dimensions, and not on dynamical characteristics 
of the system, apart from the presence of complete chaos. 
The effect of dimension on entanglement saturation has 
been studied by keeping the dimension of the first sub- 
system constant at N and increasing the dimension M 
of the second subspace from M ^ N to some large value. 
Thus we may think of the second spin as tending towards 
a complex bath with a quasi-continuous spectrum. It is 
observed that the entanglement saturation increases with 
M and finally gets saturated at the maximum possible 
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FIG. 2. The spectral average of the entanglement present 
in eigenstates of Ut {k — 9,e — 10) as a function of Q = M/N, 
where N = 2ji -1-1 = 33. Solid triangles are kicked top results 
with parity symmetry (ai = a2 = 0) and solid circles are the 
corresponding results without symmetry (qi — a2 — 0.47). 
Solid squares are the result of corresponding RMT Monte 
Carlo simulations and solid line is the theoretical curve Eq. 
(^). Horizontal line is the maximum possible entanglement 
(In(Af)). Inset shows the behaviour of the linear entropy. 




FIG. 3. Distribution of the eigenvalues of the RDMs 
of coupled kicked tops, averaged over all the eigenstates 
(TV — 2ji -I- 1 = 33). Solid curves corresponds to the the- 
oretical distribution function Eq. (^). 



kinematical limit, as shown in Fig. ^. For example the 
von Neumann entropy starting from ln(0.6A^) increases 
asymptotically to In N, while the linear entropy starting 
from I - 2/TV tends to 1 - 1/TV. 

We can develop a complete analytical understanding 
of these limits via RMT modeling. A typical stationary 
state of a quantum chaotic system shares properties of 
the eigenvectors of random matrices. Let us assume that 
some product basis has been used to write components 
o-nm of any state, which is real for stationary states of 
time reversal symmetric systems. Writing a„m as the 
rectangular N x M matrix A, the M dimensional RDM 
is A^A while the other RDM is the N dimensional AA'^ . 
The assumptions of quantum chaos, we have just seen, 
imply that A can be taken to have random independent 
entries, a member of the Laguerre ensemble. The RDMs 
then have the structure of correlation matrices , from 
where we directly use results for the density of states. 
Such matrices have also been studied since the early days 
of RMT as they have a non-negative spectrum |^ . The 
distribution of the eigenvalues of such matrices is known 
and thus this is the distribution of the eigenvalues of 
RDMs. The density of the eigenvalues of the RDM pi is 
given by 

max ^){X A„ij,i) 

= ^ A 

A=4(l + ^±;^), (7) 

where A G [Xmin, Xmax], Q = M/N and Nf{X)dX is the 
number of eigenvalues within A to A -|- dX. This has been 
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derived under the assumption that both M and N are 
large. Note that this predicts a range of eigenvalues for 
the RDMs that are of the order of 1/7V. For Q 7^ 1, 
the eigenvalues of the RDMs are bounded away from the 
origin, while for Q = 1 there is a divergence at the ori- 
gin. All of these predictions are seen to be borne out in 
numerical work with coupled tops. 

Fig. II shows how well the above formula fits the eigen- 
value distribution of actual reduced density matrices. 
This figure also shows that the probability of getting an 
eigenvalue outside the range [Xmim^max] is indeed very 
small. The sum in Sy can be replaced by an integral over 
the density /(A): 



S^ 



/(A)AlnAdA = ln(77V) 



(8) 



The integral in 7 can be evaluated to a generalized hy- 
pergeometric function and the final result is : 
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When the two subsystems are of equal dimension, that 
is Q = 1, then above expression gives 7 = exp(— 0.5) ~ 
0.6 and so the corresponding von Neumann entropy is 
ln(0.6A^). This is also the saturation value obtained in 
previous numerical work for the stationary states and 
time evolving states of a coupled standard map re- 
flecting universality. In another extreme case, when the 
Hilbert space dimension of the second subsystem is very 
large compared to that of the first, that is Q ^ 1, then 
7 ~ 1 and hence the corresponding von Neumann entropy 
is In(A^). Therefore, the analytical formulation based on 
RMT is able to explain the saturation behaviour of the 
von Neumann entropy or quantum entanglement very ac- 
curately 

Fig. g also compares the Eq. (||) to both RMT simu- 
lations and kicked top results. We expect that the devi- 
ations of the quantum calculations are due to finite size 
effects. The presence of parity symmetry results in a 
somewhat smaller entanglement, as seen in this figure, a 
fact that needs further study. For time evolving states 
and stationary states of system without time reversal 
symmetry the RDMs are complex Hermitian matrices. 
The entanglement bounds discussed here are also valid 
for these cases as the entanglement depends only on the 
density of states of the RDMs. However, spectral fluctu- 
ations of the RDMs (such as their NNSD) corresponding 
to these states can be distinct. Indeed, in the correlation 
matrix approach to atmospheric data, such a difference 



has been recently noted |10|. 

The linear entropy can also be derived as above, but 
using direct RMT results, without taking recourse to the 
distribution above, is also possible in this case. Thus we 
may write: 



N M 

Trp^ — ^ ^ ajaakaakpa-jp 

j,k=l a, [3=1 



(10) 



Substituting RMT ensemble average values [ll| of 



1/[MN{MN + 2)], {a 



3/[MN{MN + 2)] and (ajaakaCikisajis) = 0, where j ^ k 
and a 7^ /3 in the above expressions, we find that: 



1 - Tipt = 1 



M + N+1 
MN + 2 ' 
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When the dimension of the two subsystems are equal, 
that is M = N, then in the large N hmit Sr^ 1~2/N. 
This is the saturation value of the linear entropy approx- 
imately obtained in case of time evolving states of cou- 
pled kicked tops. Similarly, when the Hilbert space di- 
mension of the second subsystem is very large compared 
to the dimension of first subsystem, that is M ^ A^, 
Sr ~ 1 — This is the maximum possible value of 

linear entropy. 

Finally as an almost trivial corollary we note that the 
entanglement of formation of a time evolving density ma- 
trix is also statistically bounded to ln(77V) , as each pure 
state belonging to an ensemble representation evolves to 
this entanglement under quantum chaos. To summarize, 
we have pointed out that the eigenvalue distribution of 
reduced density matrices of composite quantum chaotic 
bipartite systems are universal, and shown that there ex- 
ists a typical value of quantum entanglement that quan- 
tum chaos engenders. This value is the maximum we 
may expect typical unentangled initial states to be able 
to reach under generic interactions. If we already had 
maximally entangled states, then chaos can disentangle 
this state to just such an extent as to coincide with this 
generic value. 
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